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STATEMENT  OF  SCOPE  AND  PURPOSE 

This  paper  presents  new  computational  methods  that  facilitate  digital  com- 
puter analysis  of  some  important  military  operations  research  problems.   Even  though 
combat  between  two  military  forces  is  a  complex  random  process,  deterministic 
Lanchester-type  differential  equations  [1]  are  commonly  used  in  defense  planning 
studies.   A  so-called  attrition-rate  coefficient  in  such  a  combat  model  represents 
the  fire  effectiveness  of  a  weapon-system  type  against  a  particular  target  type, 
i.e.  its  effective  firepower.   Time-dependent  attrition-rate  coefficients  are  used 
to  model  temporal  variations  in  firepower  on  the  battlefield.   For  such  a  variable- 
coefficient  Lanchester-type  combat  model,  we  present  a  simple  numerical  procedure 
that  allows  one  (without  having  to  explicitly  solve  the  equations)  not  only  to  predict 
battle  outcome  but  also  to  parametrically  tradeoff  quality  versus  quantity  of  two 
opposing  weapon  systems. 
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ABSTRACT 

This  paper  presents  a  simple  numerical  procedure  for  determining  the 
parity-condition  parameter  for  Lanchester-type  combat  between  two  homogeneous 
forces.   The  combat  studied  is  modelled  by  Lanchester-type  equations  of  modern 
warfare  with  time-dependent  attrition-rate  coefficients.   Previous  research  has 
shown  that  the  prediction  of  battle  outcome  (in  particular,  force  annihilation) 
without  having  to  spend  the  time  and  effort  of  computing  force-level  trajectories 
depends  on  a  single  parameter,  the  so-called  parity-condition  parameter,  which 
only  depends  on  the  attrition-rate  coefficients.   Unfortunately,  previous  research 
did  not  show  how  to  generally  determine  this  parameter.   We  present  general 
theoretical  considerations  for  its  numerical  noniterative  determination.   This 
general  theory  is  applied  to  an  important  class  of  attrition-rate  coefficients 
(offset  power  attrition-rate  coefficients) .   Our  results  allow  one  to  study 
such  variable-coefficient  combat  models  almost  as  easily  and  thoroughly  as 
Lanchester's  classic  constant-coefficient  model. 


1.   Introduction 

As  a  consequence  of  pioneering  work  by  F.  W.  Lanchester  [16]  done  about  the  time 

+ 
of  World  War  I,  military  operations  analysts  have  used  simplified  deterministic  dif- 
ferential-equation models  to  develop  insights  into  the  dynamics  of  combat  from  about 
the  end  of  World  War  II  (see,  for  example,  [1],  [7-8],  [13],  [24-28]).   The  advent  of 
the  modern  high-speed  digital  computer  has  made  feasible  the  development  and  use  of 
quite  complicated  versions  of  such  Lanchester-type  models  as  practical  defense  planning 
tools  [6].   Thus,  today  militarily  realistic  computer-based  Lanchester-type  models  of 
quite  complex  combat  systems  have  been  developed.   Such  models  currently  exist  for 
almost  the  entire  spectrum  of  combat  operations,  from  combat  between  battalion-sized 
[9]  and  division-sized  [10]  units  to  theater-level  operations  [12,14]. 

A  simple  combat  model,  however,  may  yield  a  clearer  understanding  of  important 
interrelationships  that  are  difficult  to  perceive  in  a  more  complex  model,  and  such 
insights  can  provide  valuable  guidance  for  more  detailed  computerized  investigations 

(see  [7,28]).   In  this  paper  we  present  a  new  important  numerical  procedure  that 

§ 
facilitates  parametric  analysis  of  battle  outcomes  for  such  simplified  Lanchester- 
type  models  of  combat  between  two  homogeneous  forces  with  temporal  variations  in  each 
side's  fire  effectiveness.   Previously,  such  battle-outcome  information  could  only  be 
readily  obtained  from  constant-coefficient  models."   These  results  are  not  only 
significant  in  their  own  right  but  are  also  useful  in  the  quantitative  analysis  of 
time-sequential  combat  strategies  (see,  for  example,  [19-20]). 

It  is  important  for  the  military  operations  analyst  to  have  a  clear  understand- 
ing of  how  force-level  and  weapon-system-performance  parameters  interact  to  determine 


t 
Corresponding  stochastic  formulations  are  for  all  practical  purposes  analytically 

intractable  (see  Note  1  on  p.  65  of  [22]). 

+Also  frequently  called  differential  models  of  combat  [10] . 

§ 
In  particular,  the  parametric  examination  of  force-annihilation  prediction. 

"S.  Bonder  [1-3]  has  emphasized  the  deficiencies  of  constant-coefficient  models 
(see  Section  3  below) . 


a  battle's  outcome.   Such  knowledge  is  particularly  useful  in  weapon-system  and  force- 
level  planning  activities  for  defense  planning.   S.  Bonder's  [1,3,7]  pioneering  work 
on  methodology  for  the  evaluation  of  military  systems  (particularly  mobile  systems 
such  as  tanks,  mechanized  infantry  combat  vehicles,  etc.)  provides  a  motivation  for 
interest  in  variable-coefficient,  deterministic,  Lanchester-type  combat  models  such 
as  we  consider  in  this  paper.   He  has  stressed  (see  pp.  30-31  of  [7])  the  importance 
of  analytical  solutions  to  such  models  for  developing  insights  into  the  dynamics  of 
combat  by  portraying  the  relation  between  various  factors  in  the  combat  attrition 
process  and  the  surviving  numbers  of  forces  and  for  facilitating  sensitivity  and 
other  parametric  analyses  (see  [5]).   Unfortunately,  as  work  by  Bonder  and  Farrell 
[7]  and  Taylor  [18,22]  shows,  the  analytical  (i.e.  infinite  series)  solution  to 
variable-coefficient  equations  generally  by  itself  provides  little  information 
about  battle  outcome  because  of  its  complexity.   Therefore,  one  must  seek  new  ways 
for  developing  insights. 

Taylor  and  Comstock  [24]  have  given  results  that  allow  one  to  predict  battle 
outcome  (in  particular,  force  annihilation  )  in  theory  without  having  to  spend  the  time 
and  effort  of  computing  force-level  trajectories.   To  be  computationally  practical, 
however,  their  results  require  the  determination  of  the  so-called  parity-condition 
parameter  ("the  enemy  force  equivalent  of  a  friendly  force  of  unit  strength"),  which 
depends  on  only  the  model's  attrition-rate  coefficients.   They  analytically  determine 
the  parity-condition  parameter  for  power  attrition-rate  coefficients  with  "no  offset," 
which  allow  one  to  model  combat  between  two  weapon  systems  with  the  same  maximum 


"^Especially  since  one  frequently  uses  models  that  are  so  complicated  that  trends  are 
not  directly  discernible  without  extensive  (and  costly)  computer  runs. 

t 

I.e.  without  explicitly  computing  force- level  trajectories. 

^Bonder  and  Honig  [8]  point  out,  however,  that  force  annihilation  may  not  always 
be  the  best  criterion  for  evaluating  military  operations.   See  pp.  192-242  of  Bonder 
and  Farrell  [7]  for  a  detailed  Lanchester-type  analysis  of  an  attack  scenario  for 
which  other  "end  of  battle  conditions"  play  the  principal  role.  Nevertheless,  it 
is  of  considerable  interest  (especially  for  developing  insights  into  the  dynamics 
of  combat)  to  be  able  to  easily  predict  the  occurrence  of  force  annihilation. 


effective  range  but  different  range  dependencies  for  each  system's  fire  effective- 
ness (see  also  [23]).   It  is  the  purpose  of  this  paper  to  show  how  to  determine  the 
parity-condition  parameter  in  other  cases,  in  particular  for  power  attrition-rate 
coefficients  with  "positive  offset,"  which  allow  one  to  model  such  combat  between 
weapon  systems  with  different  maximum  effective  ranges.   Our  results  allow  one  to 
study  in  general  such  variable-coefficient  combat  models  almost  as  easily  and 
thoroughly  as  Lanchester's  classic  constant-coefficient  model. 

The  organization  of  this  paper  is  as  follows.   We  first  review  Lanchester-type 
equations  of  modern  warfare,  especially  variable-coefficient  formulations.   Next  we 
review  force-annihilation-prediction  conditions  for  such  models  and  show  how  to  use 
knowledge  about  the  parity-condition  parameter  for  one  set  of  attrition-rate 
coefficients  to  numerically  determine  it  in  related  cases  of  interest.   This  general 
theory  is  then  applied  to  the  important  case  of  offset  power  attrition-rate 
coefficients,  with  detailed  numerical  examples  being  given. 

2.   Lanchester's  Classic  Combat  Formulation. 

+ 
F.  W.  Lanchester  [16]   hypothesized  in  1914  that  combat  between  two  military 

forces  could  be  modelled  by 

dx/dt  =  -  ay  ,  dy/dt  =  -  bx  ,  (1) 

with  initial  conditions 

x(t  -  0)  =  xQ  ,  y(t  =  0)  =  yQ  ,  (2) 

where   t  =  0  denotes  the  time  at  which  the  battle  begins,  x(t)   and  y(t)   denote 
the  numbers  of   X  and  Y   at  time   t,  and   a   and  b   are  nonnegative  constants  which 
are  today  called  Lanchester  attrition-rate  coefficients  and  represent  each  side's 
fire  effectiveness.   We  will  refer  to  (1)  as  Lanchester's  equations  of  modern 


+See  also  p.  45  of  [22]. 

The  equations  are  only  valid  for   x,  y  >  0.   For  example,  the  first  becomes 
dx/dt  =  0  when  x  =  0. 


warfare.^  Various  sets  of  physical  circumstances  have  been  hypothesized  to  yield 
them:   for  example,  (a)  both  sides  are  aimed  fire  and  target  acquisition  times  are 
constant  [27],  or   (b)  both  sides  use  area  fire  and  a  constant  density  defense  (see 
p.  345  of  [13]). 

From  (1)  Lanchester  deduced  his  famous  square  law 

b(x2  -  x2(t))  =  a(y^  -  y2(t))  .  (3) 

Consider  now  a  battle  terminated*  by  either  force  level  reaching  a  given  "breakpoint" 

BP  BP 

for  example,  Y  wins  when  xf  -  x(tf)  =  xfip  =  fx  xQ  but  yf  >  yBp  =  fy  yQ,  where 

denote  final  values  and  x    denotes  X's  breakpoint  which  is  a  given 


tf,  xf,  yf  denote  tinai  values  ana  xBp 
.B 

:x 


fraction  fBP  of  his  initial  strength.   It  follows  from  (3)  that 


xQ    U  -  (ff  )2} 

Y  will  win  if  and  only  if    —  <  /— RP  ?   >  ^ 

y0   ^b{l  -  iff)2} 

BP     BP 
which  for  a  f ight-to-the-f inish  (i.e.   fx  =  fy  =  0)   becomes  the  classic  result 

x0    [a 

Y  will  win  a  f  ight-to-the-f  inish  if  and  only  if   —  <  yjr   •     (5) 

y0 
Unfortunately,  no  relationship  similar  to  (3)  holds  in  general  for  variable 
attrition-rate  coefficients.5   This  paper,  nevertheless,  shows  how  (5)  generalizes 
in  these  cases."  Recalling  that  the  time  history  of  the  X  force  level  is  given  by 

x(t)  =  x-  coshv/ab  t  -  yn  >J"r  sinh/ab  t  ,  (6) 

we  see  that  the  battle  trajectories  depend  on  the  two  weapon-system-performance 

Other  forms  of  Lanchester- type  equations  appear  in  the  literature  [13,18],  but  we 
will  not  consider  these  here. 

The  modelling  of  battle  termination  is  a  problem  area  in  contemporary  defense 
planning  studies  (see  pp.  524-525  of  [25]). 

§ 
Except  when  a(t)/b(t)  =  constant  (see  p.  48  of  [22]). 

"So  far  we  have  not  been  able  to  generalize  (4). 

4 


parameters:   (I)  the  intensity  of  combat   v/ab~,  and   (II)  the  relative  fire  effectiveness 
a/b.   Only  the  relative  fire  effectiveness,  however,  determines  the  battle's  outcome 
[see  (4)  and  (5)  above]. 

3.   Variable  Attrition-Rate  Coefficients 

Bonder  [2]  has  pointed  out  that  in  many  cases  (for  example,  in  the  case  of  mobile 
weapon  systems)  the  validity  of  the  assumption  of  constant  attrition-rate  coefficients 
is  open  to  question  (see  also  [1,3,7]).   Thus,  we  consider 

dx/dt  =  -  a(t)y,  dy/dt  =  -  b(t)x  ,  (7) 

where  a(t)   and  b(t)   denote  time-dependent  attrition-rate  coefficients.   We  assume 

that  a(t)   and  b(t)   are  defined,  positive,  and  continuous  for   t„  <  t  <  +  °°  with 

t 
tn  —   0>   We  also  assume   that   a(t),  b(t)  £  L(tn,T)   for  any  finite  T.   We  further 

take  a(t)   and  b(t)   to  be  given  in  the  form  a(t)  =  k  g(t)  ,  b(t)  =  k,  h(t) ,  where 

k  ,  k,   are  positive  constants  chosen  so  that  a(t)/b(t)  =  k  /1c   when  g(t)  =  h(t). 

Analogous  to  the  constant-coefficient  case  [see  discussion  after  (6) ]  ,  we  have  the 

two  weapon-system-performance  parameters:   (I)  the  intensity  of  combat,  I(t)  = 


/a(t)  b(t);  and  (II)  the  relative  fire  effectiveness,  R(t)  =  a(t)/b(t).  We  accord- 
ingly introduce  the  combat-intensity  parameter  X  and  the  relative-fire-effective- 
ness parameter  A   defined  by 


AI  =  AaK    >       and  AR  =  ka/kb   •  (8) 

Two  significant  developments  in  the  Lanchester  theory  of  combat  during  the 
1960's  were  the  development  of  methodology  for  (a)  the  prediction  of  Lanchester 
attrition-rate  coefficients  from  weapon-system-performance  data  by  S.  Bonder  [2,4], 
and   (b)  the  (maximum  likelihood)  estimation  of  such  coefficients  from  Monte  Carlo 
simulation  output  by  G.  Clark  [11].   Both  these  developments  and  others  have 


+ 
For  convenience,  we  introduce  the  notation  that  a(t)  £  L(t„,T)   means  that 

T 
/  a(t)dt  exists  (and  is  given  by  a  finite  quantity).   From  our  assumptions  about 

t0  T 

a(t)   and  b(t),  it  follows  that  a(t)  t   L(tQ,T)  implies  that  J*t  a(t)dt  =  +  °°  . 

See  [22]  for  further  references. 


generated  interest  in  the  model  (7)  and  facilitated  its  application  (and  that  of 
its  generalization  to  combat  between  heterogeneous  forces  [7])  to  defense  planning 
studies. 

A  large  class  of  tactical  situations  of  interest  can  be  modelled  with  the 
following  general  power  attrition- rate  coefficients  [7,22,24] 

a(t)  =  k  (t  +  C)y,        and       b(t)  =  k^t  +  C  +  A)v,         (9) 
where  A,  C  >^  0.   We  will  call  A  the  offset  parameter,  since  it  allows  us  to  model 
(with  y,  V  _>  0)  battles  between  weapon  systems  with  different  maximum  effective 
ranges.   We  will  call  C   the  starting  parameter,  since  it  allows  us  to  model  (again 
with  y,  v  >_  0)  battles  that  begin  within  the  minimum  of  the  maximum  effective  ranges 
of  the  two  systems.   For  example,  let  us  consider  Bonder's  [1,3]  model  of  a  constant- 
speed  attack  on  a  static  defensive  position  (see  also  [18,22]).   Then  we  have 

dx/dt  =  -  a(r)y  ,  dy/dt  =  -  3(r)x  ,  (10) 

where  r(t)  =  R~  -  vt  denotes  the  distance  (range)  between  the  two  opposing  forces, 
R-.  denotes  the  battle's  opening  range,  v  >  0  denotes  the  constant  attack  speed, 

for  r>Ra, 


(11) 
aQ(l-r/Ra)y         for  0  <  r  <  Ra> 

y  J>  0,  and  R   denotes  the  maximum  effective  range  of  Y's  weapon  system.   Similarly 
for  $(r),  with  exponent  v  >^  0.   In  (11)  the  parameter  y  allows  us  to  model  the 
range  dependence  of  Y's  fire  effectiveness  (see  Fig.  1).   The  offset  and  starting 
parameters  are  given  by 

A  =  (Rg  -  Ra)/v  ,        and        C  =  (Ra  -  RQ)/v  ,  (12) 

and  the  assumption  A,  C  _>  0  implies  that  RR  >.  ^  ^.  Rn«   From  considering  (12) 
and  Fig.  2,  the  reader  should  have  no  trouble  understanding  our  terminology  for  A 
and   C. 

The  time  history  of  the  X  force  level,  i.e.  the  solution  x(t)   to  (7),  is 
given  by  [22] 
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x(t)  =  xQ{CY(0)Cx(t)-SY(0)Sx(t)}  -  yQAR  {Cx(0)Sx(t)-Sx(0)Cx(t)},        (13) 

where  the  hyperbolic-like  general  Lanchester  functions  (GLF)   Cy(t)   and  Sy(t)   are 
linearly  independent  solutions  to  the   X  force-level  equation 


dt2 


-{^yi}f  - -<t>  «.<t)x  -  o ,  (H) 


with  initial  conditions 


cx(t  =  tQ)  =  1  ,  sx(t  =  tQ)  =  0  , 

(15) 
{[l/a(t)]dCx/dt(t)}t=t  =  0  ,       {[l/a(t)]dSx/dt(t)}t=t  =  1/A^  , 

where   t„   denotes  the  largest  finite  time  at  which  a(t)   or  b(t)   ceases  to  be 
defined,  positive,  or  continuous.   For  example,  t_  =  -C  for  the  general  power  attti- 
tion-rate  coefficients  (9) .   The  time  history  of  the  Y   force  level  may  be  similarly 
obtained,  with  C  (t)   and  S  (t)   being  analogously  defined  for  the  corresponding 
Y   force-level  equation. 

For  the  numerical  determination  of  the  parity-condition  parameter,  it  is  con- 
venient to  introduce  a  new  independent  variable  s  defined  by 

t 
s  =  KAI  ft     g(a)  da  ,  (16) 

where  the  parameter  K  is  to  be  chosen  to  simplify  the  form  of  J(s)   given  by  (18) . 
We  denote  s(t  =  0)   as   s  ,  and  then  sn  >^  0  if  and  only  if   t_  £  0.   The  sub- 
stitution (16)  transforms  (14)  into  the  normal  form  [17] 

d2x 

-  J(s)x  =  0  ,  (17) 


ds2 


where 


'<«■-;*{!£}   • 


K 
and  t  =  t(s).   We  also  define  the  normal-form  hyperbolic-like  GLF  c Y(s)   and 


X 


s  (s),  which  satisfy  (17)  and  the  initial  conditions 


cx(s   =  0)    =  1    ,  sx(s  -  0)    =   0    , 


dcx/ds(s  =  0)   =  0    ,  dsx/ds(s  =  0)   =  1    . 


(19) 


It  follows  that 

cx(s(t))  =  Cx(t)  ,    and       sx(s(t))  =  KSx(t)  .  (20) 

4.   Force-Annihilation-Prediction  Conditions 

THEOREM  1  (Taylor  and  Comstock  [24]):   Assume  that  either  a(t)  t   L(0,+  °°)   or 

b(t)  t   L(0,+  °°) .   Then  the  X  force  will  be  annihilated  in  finite  time  if  and  only  if 


X  JT    1  CX(0)  "  Q*SX(Q) 


>\   -t ~  1  .  (2D 

yo       (q  Cy(0)-SY(0) 

* 


where  the  parity  condition  parameter  Q   i£  unique  and  given  by 

sx(t)   ^ 

lim  C   (t~)   ~     *    '  (22^ 

f*-+«>  uXkt;    Q 

Remark  1:  We  also  have   lim  {S  (t)/C  (t)}  =  Q 

Remark  2:  When  t  =  0,  (21)  simplifies:   X  will  be  annihilated  in  finite  time  if 


and  only  if  xQ/y0  <  »^/Q   • 


Remark  3:   The  result  (22)  suggests  a  numerical  procedure  for  approximately  deter- 
mining  the  parity-condition  parameter  Q  :  we  may  approximate  the  parity-condition 
parameter  Q   by  Q  =  l/{SY(t) /Cv(t) } ,  where   t  is  a  "suitably  large"  value  of   t. 

* 

In  other  words,  we  may  estimate  Q   simply  by  picking  a  large  value  for   t   (we 
denote  this  selected  large  value  by  t),  computing  Sv(t)   and  Cv(t) ,  and  then  form- 

A  A 

ing  their  ratio.   Our  estimate  for  Q   is  then  given  by  Q  =  1/{S  (t)/Cv(t)}.   The 

A  A 

only  problem  is  that  we  don't  know  how  large  to  take   t  for  "satisfactory"  estima- 

•k  *  -k  a  a 

tion  of  Q  :   There  is  an  estimation  error,  E(t)  =  Q  -  Q(t),  which  depends  mono- 

tonically  on  t,  and  a  priori  we  don't  know  how  large  this  error  is.   The  present 

paper  develops  a  bound  on  the  magnitude  of  this  error,  and  our  new  error  estimate 
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allows  the  goodness  of  approximation  to  be  easily  evaluated  in  many  cases  of  interest. 
We  may  also  determine  the  parity-condition  parameter  with  the  normal-form  hyper- 
bolic-like  GLF,  since  lim    ^{s  (s) /c  (s)  }  =  1/Z   =  K/Q  ,  where   Z    is  called  the 
modified  parity-condition  parameter.   In  fact,  we  will  find  it  more  convenient  to  do 
so.   With  this  in  mind,  let  us  introduce  the  Y-functions   c  (s)   and  s  (s)  [corres- 
ponding to  cy(s)   and  s  (s)]   defined  by 

dc^/ds  =  J(s)sx  ,  dsy/ds  =  J(s)cx  ,  (23) 

with  initial  conditions 

(^(s  =  0)  =  1  ,  sy(s  =  0)  =  0  .  (24) 

It  follows  that   c  (s)   and  s  (s)   are  linearly  independent  solutions  to  the  modified 


Y  equation 


and 


_d_ 
ds 


iikfs}-y-°'  <»> 


cy(s(t))  =  Cy(t),  sY(s(t))  =  (1/K)  Sy(t)  .  (26) 

In  terms  of  the  new  time  variable  s  defined  by  (16) ,  Theorem  1  reads  as  follows 

THEOREM  2:   Assume  that  either  a(t)  t   L(0,  +  °°)   or  b(t)  k   L(0,  +  °°)  .   Then  the  X 
force  will  be  annihilated  in  finite  time  if  and  only  if 


^0   ^R    \   W  ■  Z  SX(S0> 

yo  <  K  \  z\v( 


(27) 


rvS0)     SY(SQ) 


where  the  modified  time  variable  s   is  given  by  (16),  and  c„(s) ,  s  (s) ,  cy(s) ,  and 
sy(s)   denote  the  normal- form  hyperbolic- like  GLF.   The  modified  parity-condition 

parameter   Z   _is_  unique  and  given  by 

s  (s) 
lim   -*—  =  -\  .  (28) 

s  -»•+  »   CX^S;    Z 


We  observe  that  .      . 


Q  =  KZ   ,  (29) 


and  lim    ,  isv(s) /cv(s)  }  =  Z  .   When  (27)  holds,  the  time  to  annihilate  X,  denoted 

S  ->  +00   I        I 

X  X 

as   t  ,  is  determined  by  x(t  )  =  0.   If  we  denote  the  quotient  of  the  two  normal- 

3.  3- 

form  hyperbolic-like  GLF  s  (x)   and  c  (s)  as  n„(s) ,  then  it  follows  from  (13)  that 

X         x        X 

,  ,JC„    {X0CY(S0)  +yQ(^R/K)SX(s0)> 

nx(s(ta))  = ,  (30) 

{xQsY(s0)  +  y0(A^/K)cx(s0)} 


where 


nx(s)  =  sx(s)/cx(s)  .  (31) 


5.   Determination  of  the  Parity-Condition  Parameter. 

We  will  now  show  how  knowledge  about  the  modified  parity-condition  parameter  Z 
for  one  pair  of  attrition-rate  coefficients,  a(t)  and  b. (t) ,  allows  us  to  determine 
Z   for  a  related  pair,  a(t)   and  b(t).   With  this  in  mind,  let  us  denote  cy(s) 
corresponding  to  a(t)   and  b(t)   as  cy(s;a,b),  and  similarly  for  s   and  r\    . 
In  other  words,  we  will  now  write  (31)  corresponding  to  the  attrition-rate  coefficients 
a(t)   and  b(t)   as 

r)x(s;a,b)  =  sx(s;a,b)/cx(s;a,b)  .  (32) 

f 
In   this   notation,  we  will  write    (28)    as 

lim       r^(s;a,b)   =  1/Z*[a,b]    .  (33) 

s  -*■  +°° 

Our  main  result  is  Theorem  5,  which  gives  an  error  estimate  for  the  approximation 
that  we  propose  for  Z  .   The  theoretical  basis  for  Theorem  5  is  given  by  Theorem  4, 
which  (in  turn)  is  a  consequence  of  Theorem  3.   The  proof  of  Theorem  3  follows  along 
the  lines  of  well-known  arguments  (see  p.  225  of  [15]). 


t  * 

We  use  the  notation  Z  [a,b]   to  show  that  the  modified  parity-condition  parameter  is 

a  functional  (i.e.  a  function  for  which  the  independent  variables  themselves  are 

functions),  which  depends  on  only  the  attrition-rate  coefficients  a(t)   and  b(t). 

In  other  words,  the  attrition-rate  coefficients  are  functions  defined  for  t~  _<  t  < +°°  , 

and  the  parity-condition  parameter  depends  on  these  entire  functions  (and  not 

merely  particular  values  of  them) . 
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THEOREM  3  (Comparison  Theorem):   Let   x(t)   and   x  (t)   satisfy 

IFUof  }-"(t)x-o,       A{^^lj.bi(t)Xi  =  0, 

with  initial  conditions 

x(t  =  tQ)  =  a  ,  x±(t  =   tQ)  =  a  , 

{[l/a(t)]dx/dt(t)}    =  3  ,  {[l/a(t)]dx  /dt(t)}    =  6  , 

0  0 

where  a(t)  >  0  and  b  (t)  <  b(t)   for  all   t  >  tQ.   Then  x  (t)  <  x(t)   for  all 

t  >  tn  as  long  as  x(t)  >  0. 

The  basic  theoretical  result  upon  which  our  numerical  determination  of   Z 
is  based  is 
THEOREM  4.   Assume  that   b  (t)  <  b(t)   for  all   t  >  t„.   Then 

nx(s;a,b)  <  1/Z*[a,b]  <  nx(s;a,b)  +  {  (1/Z*[a,b1] )  -r^s  ;a  ,1^)  }  .      (34) 
PROOF :   We  observe  that  [24]   r|v(s;a,b)   satisfies  the  differential  equation 

drix/ds(s;a,b)  =  l/{cx(s  ;a,b)  }2  ,  (35) 

with  r)v(s  =  0;a,b)  =  0,  and  similarly  for  nv(s;a,b1).   Theorem  3  (the  comparison 
theorem)  yields  that  c  (s;a,b)  >  cv(s;a,b1)   for  all  s  >  0.   Thus,  for  all  s  >  0 

X  X        J_ 

dnx/ds(s;a,b)  <  dnx/ds(s ;a ,b1>  , 
whence  integration  between  0  and  s  yields  the  desired  result.  Q.E .D. 

Similar  to  the  observations  made  in  Remark  3  above,  we  observe  that  (33) 
suggests  that  we  estimate   Z  [a,b]   with   Z  defined  by 

Z(s;a,b)  =  l/nx(s;a,b)  ,  (36) 

where  s  denotes  a  suitably  chosen  value  for  s.   Moreover,  from  (35)  we  see  that 
n.„(s;a,b)   is  a  strictly  increasing  function  of   s  so  that  the  larger  we  take  s 
in  (36),  the  better  our  approximation  becomes.   The  only  problem  (see  Remark  3)  is 
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that  a  priori  we  don't  know  how  large  to  take  s   for  "satisfactory"  estimation 

* 
of  Z  .   Theorem  5,  however,  tells  us  exactly  how  large  to  take  s. 

THEOREM  5  (Error  Estimate  for  Approximation):   Assume  that  b..  (t)  <  b(t)   for  all 

*  * 

t  >  tn.   Let  f^Cs)   denote  the  fractional  error  made  in  the  estimation  of  Z  [a,b] 

(J    hi — - — - — 

by  Z ( s ; a , b ) ,  i.e. 

f_ (s)  =  Z(s;a,b)  -  Z*[a,b]  (37) 

Z  [a,b] 
Then 

0  <  fE(s)  <  {(1/Z*[a,b1])  -  nx(s;a,b1)}  ■  Z(s;a,b)  .  (38) 

PROOF :   The  theorem  follows  by  simple  algebraic  manipulation  after  setting  s  =  s 
in  (3A)  and  using  (37).  Q.E.D. 

Thus,  we  have  presented  a  method  for  numerically  determining  Z  [a,b].   We 
simply  pick  a  large  value  for  s   (we  denote  the  selected  value  as  s) ,  compute 
s  (s)   and  cv(s) ,  and  then  compute  the  estimate  Z(s;a,b)   according  to  (36). 

A.  A. 

Theorem  5  allows  us  to  know  the  accuracy  of  our  approximation,  which  can  be  improved 

* 
by  taking  s   larger.   Thus,  we  can  numerically  determine  Z  [a,b]   to  any  specified 

degree  of  accuracy  once  Z  [a,b- ]   is  known.   In  the  next  section  we  apply  this 

theory  to  the  analysis  of  battles  modelled  with  offset  power  attrition-rate 

coefficients. 


6.   Application  of  Theory  to  Offset  Power  Attrition-Rate  Coefficients 

In  the  application  of  Theorems  4  and  5,  two  pairs  of  attrition-rate  coefficients 
are  involved:   one  pair  for  which  the  modified  parity-condition  parameter  is  known 
[denoted  as  a(t)   and  b1(t)],  and  one  for  which  it  is  to  be  determined  [denoted 
as  a(t)   and  b(t)].   Accordingly,  we  rewrite  (9)  with  A  >  0  as 

a(t)  =  k  (t  +  C)y,  and         b(t)  =  k,  (t  +  C  +  A)V,     (39) 
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where  (as  before)  C  >_  0.   We  will  refer  to  these  coefficients  (39)  for  which  A  >  0 
as  power  attrition- rate  coefficients  with  "positive  offset."   If  we  choose 

K  =    [Xj/Cy+l)]21*'1    ,  (40) 

it  follows  from  (16)  that  the  modified  time  variable  s   is  given  by 

s(t)  =  [AI/(y+l)]2p(t  +  C)y+1  ,  (41) 

and  the  invariant  J(s)   of  the  normal  form  (17)  simplifies  to 

J(s;a,b)  =  J(s;Y,]i,v)  =  s6(l  +  ^)  ,  (42) 

s 

where  p  -  (y+l)/l,  a  =  l/(y+l),  3  =  (v-y)/(y+l),  y  =  A- [A  /(y+1) ]2/Z,  and  Z  =  y+v+  2. 
Here  we  have  denoted  the  invariant  corresponding  to  the  attrition-rate  coefficients 
a(t)   and  b(t)   as  J(s;y,y, v) ,  since  we  may  take  y,  y,  and  v  as  a  basis  for 
generating  the  four  parameters  a;  3S  Y,  and  v   that  explicitly  appear  in  the  right- 
hand  side  of  (42).   Furthermore,  we  will  denote  the  normal-form  hyperbolic-like  GLF 
that  correspond  to  J(s;y,y,v)   as  c  (s;y,y,v)   and  s  (s;y,y,v). 

The  known  results  [24]  that  we  use  in  the  Theorems  4  and  5  are  for  the  case 
of  power  attrition- rate  coefficients  with  no  offset  [i.e.  set  A=0  in  (9)] 

a(t)  -  k  (t  +  C)y,      and      b^t)  =  k^t  +  C)V,  (43) 

where  C  >^  0.   We  observe  that  b  (t)  <  b(t)   for  all   t  >  -C.   It  follows  that 
J(s;a,b1)  =  s^  and  [24] 


Z*[a,b1]  =  p(2p  X)  r(l-p)/r(p)  .  (44) 


Thus,  for  the  bound  on  Z  [a,b]  =  Z  (y,y,v)   given  in  Theorem  4  and  the  error  estimate 
for  our  approximation  (36)  given  in  Theorem  5,  we  have  [23] 

nx(s;a,b1)  =  p(1_2p)  T  (S)  ,  (45) 

where  S  =  2p  s    P  ,  q  =  1-p,  and  T   denotes  a  Lanchester-Clif ford-Schaf li 

f 
(LCS)  function,  which  is  analogous  to  the  hyperbolic  tangent  (see  Table  I) . 


f 
These  functions  were  introduced  in  [22]  and  redefined  for  reasons  of  force-annihi- 
lation prediction  in  [23]. 
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TABLE  I. 
LANCHESTER-CLIFFORD-SCHAFLI  FUNCTIONS 


°°         (   /:n2k 

F«(x)  =  r(0>  Jo  fc!  *<*  *  «> 

,     2(k-K») 

Ha(x)  =  r(«)  Jq  k.  $2>  „  +  iy 


T    (x)    =  H,       (x)/F    (x) 
or  1-or         a 


Relation  to  normal-form  GLF. 


cx(s)    =  Fq(S)  sx(s)    =  p(1"2p)Hp(S) 

cy(s)    =  F   (S)  sY(s)    =  p(2p"1)Hq(5) 


where     q  =   1-p     and 


S(s)    *  2pa1/(2P) 


NOTE:   For  ]i   =  V,  we  have 


(I)   cx(s)  =  cy(s)  =  Fh/2(s)  =  cosh  s, 

(II)   sx(s)  =  sy(s)  -  H^Cs)  -  sinh  s, 

(III)   nx(s)  =  n  (s)  m   T«  /2(s)  m   tanh  s. 


In  fact, 


n(s;a,b..)  =  tanh  s,        when  y  =  v  .  (46) 


We  have  thus  shown  that  the  following  theorem  holds. 
THEOREM  6 :   Assume  that  either   y  >  -1   or  v  >  -1 .   Then  for  a  battle  modelled  with 
the  offset  power  attrition-rate  coefficients  ( 39 ) ,  bounds  on  the  modified  parity- 
condition  parameter  Z  (y,y,v)   are  given  for  y  >  0  by 

nx(s;y,y,v)  <  -j— ^ <  nx(s;y ,y,v)  +  pq_p  {  -j^-  -  T  (S)  \     ,        (47) 

z  (y,y,v)  ^q;    q 

where  q  =  1-p,  S  =  2ps    p  ,  and  n(s;y,y,v)  denotes  the  quotient  of  two  normal- 

A  "' 

form  hyperbolic-like  GLF  for  the  attrition-rate  coefficients  (39),  i.e.   nv(s;y,y,v) 
=  sx(s;y,y,v)/cx(s;y,y,v). 

It  follows  from  Theorem  6  (or,  equivalently ,  Theorem  4)  that  if  we  approximate 
Z  (Y.y.v)   with  Z(s;y,y,v)   defined  by 

Z(s;y,y,v)  =  l/nx(s;y ,y ,v)  ,  (48) 

then  bounds  on  the  fractional  error  made  in  this  estimate  are  given  by 

0  <  fE(s)  <  Pq_P{r{^}-  Tq(S)|  nx(s;y,y,v),  (49) 

where   f^s)   denotes  the  fractional  error  and  is  defined  by  (37)  . 
E 

The  right-hand  inequality  in  (49)  [equivalently,  (47)]  tells  us  exactly  how 
large  to  take  s   for  the  estimation  of   Z  (y  >  0,y,v)   by  Z(s;y,y,v)   to  any 
specified  degree  of  accuracy.   The  LCS  function  T   is  involved  in  the  bound  on 
the  fractional  error  fT7(s)   in  this  estimate  when  y  ^  v.    Thus,  the  LCS  functions 
as  redefined  by  Taylor  and  Brown  [23]  yield  valuable  information  about  battles 
modelled  with  not  only  the  power  attrition-rate  coefficients  with  no  offset  (43)  but 
also  the  offset  power  attrition-rate  coefficients  (39).   Availability  of  tabulations 
of  these  LCS  functions  is  discussed  in  [23]. 


4- 

This  result  is  one  of  our  reasons  for  introducing  the  normal  form  (7). 

*As  (46)  and  Table  I  show,  T  (S)  =  tanh  s  when  y  -  V. 
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7.   Numerical  Results 

In  this  section  we  will  examine  a  couple  of  numerical  examples  to  show  how 
the  modified  parity-condition  parameter  Z   may  be  numerically  determined  and  to  show 
some  important  insights  into  the  dynamics  of  combat  that  may  be  consequently  obtained. 
In  order  to  numerically  determine  the  modified  parity-condition  parameter  for  the  off- 
set power  attrition-rate  coefficients  (39),  we  must  use  knowledge  about  how  quickly 
the  limiting  value  (i.e.   Z  [a,b  ])  of  a  hyperbolic-tangent-like  function  of  a  related 
pair  of  power  attrition-rate  coefficients  with  "no  offset"  (43)  is  reached  as  its 
argument  is  increased  [recall  Theorem  6  and  (49)].   In  Fig.  3  we  see  that  this  limit- 

ing  value,  denoted  as  Z  (y,v)  =  Z  [a,b..],  is  quite  quickly  reached:   if  one  takes 

* 
s  =  10.0,  then  Z  (y,v)   is  approximated  to  better  than  six  decimal  places  by 

Z(s;y,v)  =  l/nx(s;y,v),  where  T)        is  given  by  (45).   Experimental  computing  for 

various  values  of  y  and  v  and  comparison  with  the  known  value  (44)  for  Z  (y,v) 

bears  out  this  degree  of  accuracy  [i.e.  speed  of  convergence  of  Z(s;y,v)   to  Z  ] 

•\  /\ 
for  essentially  all  allowable  values  of  y  and  v.   Thus,  Z(s;y,v)   for  the  coef- 

ficients  (43)  has  essentially  converged  to  Z  (y,v)  when  s  =  10.0,  and  by  Theorem  6 

or  (49)  we  know  that  the  same  is  true  for  Z(s;y,y,v)   for  the  coefficients  (39). 

We  have  accordingly  generated  by  this  procedure  the  results  shown  in  Fig.  4. 
For  computing  r\     =   sy/c„,  we  have  used  the  series  solutions  shown  in  Tables  II  and 
III.   [In  Tables  II  and  III  we  have  for  convenience  denoted,  for  example,  sv(s;Y,y,v) 
simply  as  s  (s;y,v),  i.e.   sy(s;y,v)   denotes  s   corresponding  to  the  general 
power  attrition-rate  coefficients  (9)  with  exponents  y  and  v.]   The  series  were 
obtained  by  solving  (17)  by  the  method  of  successive  approximations  (see  [18]).   We 
used  these  series  instead  of  developing  approximate  solutions  by  finite-difference 
methods  because  we  did  not  have  any  error  bounds  for  the  latter. 

Let  us  now  give  an  intuitive  interpretation  of  the  curves  shown  in  Fig.  4 

of  the  modified  parity-condition  parameter  Z   plotted  versus  the  modified  offset 

parameter  Y.   In  Taylor  and  Comstock  [24]  it  is  shown  that  Z   may  be  considered 
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TABLE  II 
NORMAL-FORM  OFFSET  LINEAR  LANCHESTER  FUNCTIONS 

00  2k        k         .    /         Vj 

cx(s;l,l)    =      I     ^yr     I      Aj  U- 


2k+l     k  .  /  \j 

sx(s;1,l)   -  Jq  ^^fejyr 


j=o     M/r 


where  the  offset  coefficients  are  given  by 


A  =  1,    and  for  k  >  1 


4k(4k-2)     {J  Aj-1 


Tc   (4k-j)(4k-2-j)   1Ak-l  +  ^-l/       t0r  U  -  J  -  k 


B°  =  1,    and  for  k  >  1 


I 


k   (4k-j) (4k+2-j)    k-1    k-1  -J- 


We  have  adopted  the  convention  that  AJ  ,  B^  =  0   for  j  <  0  or  j  >  k 

k    k 


TABLE  III 
OFFSET  POWER  LANCHESTER  FUNCTIONS  FOR  y  =  1  AND  v  =  2 


m      ,2      5/4, 2k    01 

3   v      ( 5  S    }      ?k    i  /  Y  \J 
.x(s;l,2)  =  r(|)  £-A__   ^  ^£j 


k=0  k!  T(k 


00 

2,1/5  r,2, 

s 


(2  5/4,2(1*2/5)   2k 


(s;l,2)  -  (|)1/5  r(|)   I  ^ I  aJ 

k=0   k!  r(k  +  I)    J=0 


f 
where  the  offset  coefficients  are  given  by 


A  =  1,    and  for  k  >  1 


aJ      5k(5k-2)     c  .3~2    .    OAj-l    .j   -,     r    „    ,    „, 

K  =  (5k-j)(5k-2-j)  {Ak-i +  2Ak-i +  ^-iJ     for   °  <  J  <  2k 


B     =  1,  and   for     k  >  1 


*l  =    FiiT7w^J9    ^      {BJJ  +  2BJJ+BJ.}        for     0  <  j   <   2k 
,  k        (5k-j) (5k+2-j)  k-1  k-1  k-1  —       — 


We  have  adopted  the  convention  that  A~,  B^  =  0   for  j  <  0  or  j  >  2k. 


t 
to  be  the  initial  Y  force  level  that  leads  to  a  draw   (i.e.  parity  between  the 

forces)  in  the  battle  against  an  X  force  of  "unit  strength" 

dx/ds  =  -y  with   x(s  =  0)  =  1  , 

*  .  (50) 

dy/ds  =  -J(s)x         with  y(s  =  0)  =  Z   , 

where  J(s)   denotes  the  invariant  of  the  normal  form  (17).   Thus,  we  may  consider 

Z   to  be  "the  Y  equivalent  of  an  X  force  of  unit  strength"  for  the  modified 

battle  (50) .   Now  let  us  consider  the  general  power  attrition-rate  coefficients  (9) 

with  exponents  y  and  v.   As  we  did  in  Tables  II  and  III,  we  will  denote  the  corre- 

sponding  J  as  J(s;y,v)   and  Z   as  Z  (y,v)  to  stress  the  dependence  on  y  and 

V   (but  suppressing  that  on  y)  .   We  then  have  from  (42)  that  J(s;l,l)  =  1  +  y//s 

and  J(s;l,2)  =  /s    (1  +  y//s)2.   From  (44)  we  find  that  Z*(l,l)  =  1.000  and  Z*(l,2) 

=  0.806  for  y  =  0.   Observing  that  for  y  >  1  we  have  J(s;l,l)  <  J(s;l,2)   for  all 

s  ^  0,  it  is  intuitively  clear  from  (50)  and  the  interpretation  of  Z   as  a  force 

equivalent  that  we  must  have  Z  (1,1)  <  Z  (1,2)   for  all  Y  >  1  because  X  always 

has  greater  fire  effectiveness  against  Y  when  y  =  1  and  v  =  2   than  when  y  =  1 

and  V  =  1.   However,  for  y  near  zero,  the  situation  is  reversed  and  Z  (1,2)   must 

* 
lie  below  Z  (1,1)   for  y  near  zero.   Thus,  we  have  given  an  intuitive  explanation 

*  * 

of  why  Z  (1,2)   lies  below  Z  (1,1)   for  y  near  zero  but  above  it  for  y  >  1  as 

Fig.  4  shows. 

Next,  we  will  consider  numerical  results  for  a  particular  battle  to  show  some 
of  the  important  insights  that  may  be  gained  into  the  dynamics  of  combat  from  our  new 
results.   As  in  [18,22,23]  we  consider  S.  Bonder's  [1,3]  model  (10)  for  the  constant- 
speed  attack  of  mobile  forces  against  a  static  defensive  position.   We  will  focus  on 
the  new  results  given  in  this  paper  (in  particular,  the  prediction  of  battle  outcome 
from  initial  conditions  without  explicitly  computing  the  force-level  trajectories). 
Input  data  and  computed  parameter  values  are  shown  in  Table  IV.   We  will  now  consider 
two  cases:   (I)  R  =  1500  meters,  and   (II)  R  =  1250  meters. 

In  other  words,  x(s)   and  y(s)  >  0   for  all  s  e  [0,4«°)   but   lim     x(s)  =  0 
=  lim    .  y(s). 
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TABLE  IV 
PARTICULARS  FOR  THE  NUMERICAL  EXAMPLES 


1.   Input  Data 
y  =  V  =  1 


a  =  0.06  X  casualties/minute/Y  unit 


(3Q=  0.6  Y   casualties/minute/X  unit 


R  =  1500  meters,      R0=  2000  meters 
v  =  5  miles/hour 


Parameter  Values 

-3 
k  =  5.364  x  10   X  casualties/minute/Y  unit 
a 

_3 

k=  4.023  x  10   Y  casualties/minute/X  unit 

P  =  q  =  1/2 

A  =  3.728  minutes,    Y  =  0.320  (casualties 'minutes) 


1/2 


When  RQ  =  1500  meters,  we  have   C  =  0   and   s   =  0.   The  maximum  time  that 

the  battle  can  last  is   t    =  11.18  minutes,  since  at  this  time  the  advancing 

max  ° 

>'c  -k 

attackers  overrun  the  defensive  position.   In  this  case  Z  (y,u,v)  =  Z  (0.32,1,1) 

=  1.381,  so  that  Theorem  2  tells  us  that   X   can  be  annihilated   <==>  x  /y  <  0.264. 

By  (30)  the  X-force  annihilation  time  is  given  by  riv(s(t  ))  =  2.  739xn/yri .   For 

a    a  u  u 

X  + 

xn  =  10  and  yn  =  50,  we  have  nv(s  )  =  0.54772  so  that  by  the  techniques   introduced 

U  U  A   3, 

in  [23]  we  find   sX  =  0.771.   Hence,  (36)  yields   tX  =  10.25  minutes  and   rX  =  125.7 
a  a  a 

meters.   Further  results  are  given  in  Table  V. 

When  RQ  =  1250  (see  Fig.  5  of  [22]),  we  have  C  =  1.864  minutes,  s  =  0.0255 

and   t    =9.32  minutes.   In  this  case  X  can  be  annihilated  <==>  x^/y^  <  0.281 
max  0  J0 

with  the  X-force  annihilation  time  given  by  nv(sX)  =  (1.001un  +  0.009)/(0.127u  +  0.366), 

X   3.  (J  (J 

where  Uf)  =  x_/yn.   Numerical  results  are  given  in  Table  VI.   Finally,  these  parametric 
results  should  be  contrasted  to  those  previously  possible  (e.g.  compare  them  with, 
for  example,  the  single  force-level  trajectory  for  Rq  =  2000  meters  shown  in  Fig.  5 
of  [22]). 


8.   Discussion 

S.  Bonder  [1-3]  has  emphasized  the  shortcomings  of  constant-coefficient 
Lanchester-type  combat  models.   Work  by  Bonder  [1,2],  Clark  [11],  and  others  [7]  on 
the  prediction  of  Lanchester  attrition-rate  coefficients  has  generated  interest  in 
variable-coefficient  models.   Moreover,  there  is  not  only  intrinsic  interest  (see 
[1,3])  in  the  model  (7)  but  also  interest  for  obtaining  insights  into  the  behavior 
of  complex  Lanchester-type  system  models§  that  have  been  enriched  in  military  detail 


t 
These  computations  involve  the  generation  of  a  table  of   s  ,  c  ,  and  ri  for 

Y  =  0.32,  y  =  v  =  1  (see  [23]). 

See  Taylor  and  Brown  [22]  for  further  discussion  and  references. 
§For  example,  the  Bonder-IUA  model  (see  [7-9]). 
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TABLE  V 
ANNIHILATION  OF  THE  X  FORCE  AS  A  FUNCTION 
OF  THE  INITIAL  FORCE  RATIO  FOR  R  =  1500  METERS 


(x0/yQ) 

t    (minutes) 
a 

r    (meters) 
a 

0.250 

14.09 

10.25 

8.80 

t 

0.200 
0.167 

125.7 
319.4 

t 


t    =11.18  minutes  and  x£  =  x(r  =  0)  =  2.48 
max  f 


TABLE  VI 

ANNIHILATION  OF  THE  X  FORCE  AS  A  FUNCTION 

OF  THE  INITIAL  FORCE  RATIO  FOR  R  -  1250  METERS 

0 


(xQ/y0) 

a 

(minutes) 

r    (meters) 
a 

0.250 

10.87 

t 

0.200 

8.17 

154.4 

0.167 

6.93 

320.4 

t 


t    =9.32  minutes  and  x£  =  x(r  =  0)  =  1.74 
max  f 


see  [7-10,12,14]).   The  attrition-rate  coefficients  in  (7)  represent  the  fire 
effectiveness  of  the  combatants  and  allow  us  to  model  temporal  variations  in  fire 
effectiveness  on  the  battlefield.   Interest  in  the  general  power  attrition-rate 
coefficients  (9)  is  provided  by  S.  Bonder's  [1,3,8]  constant-speed  attack  modelT 
(10)-(11)  and  his  examination  of  the  range  dependence  of  attrition-rate  coefficients 
for  various  weapon  systems  (see  pp.  196-200  of  [7]). 

We  have  given  results  that  allow  one  to  study  the  variable-coefficient  model 
(7)  [especially  with  the  general  power  attrition-rate  coefficients  (9)]  almost  as 
easily  and  thoroughly  as  Lanchester's  classic  constant-coefficient  model  (1).   Taylor 
and  Comstock  [24]  (see  Theorems  1  and  2  above)  have  shown  how  to  predict  force 
annihilation  without  having  to  spend  the  time  and  effort  of  explicitly  computing 
force-level  trajectories.   Using  their  theoretical  results,  we  gave  results  in  a 
previous  paper  [23]  that  made  combat  modelled  by  power  attrition-rate  coefficients 
with  no  offset   [i.e.   A  =  0   in  (9)]  almost  as  easy  to  analyze  as  the  constant- 
coefficient  case.   The  results  of  the  paper  at  hand  allow  one  to  analyzed  combat 

§ 
modelled  by  power  attrition-rate  coefficients  with  positive  offset   [i.e.   A  >  0 

in  (9)]  just  as  conveniently. 

Theorem  1  (see  also  Theorem  2)  is  the  generalization  of  the  classic  constant- 
coefficient  result  (5)  to  cases  of  time-dependent  attrition-rate  coefficients.   How- 
ever,  one  needs  to  know  the  value  of  the  so-called  parity-condition  parameter  Q 
in  order  to  predict  force  annihilation  in  specific  instances.   In  this  paper  we  have 
presented  theoretical  considerations  (see  Section  5  above)  for  the  noniterative 
numerical  determination  of  the  parity-condition  parameter.   We  applied  our  general 
theory  to  the  specific  case  of  general  power  attrition-rate  coefficients  (9)  (see 


Thus,  the  range  between  firer  and  target  changes  during  the  engagement. 

Modelling,  for  example,  combat  between  two  weapon  systems  with  the  same  maximum 
effective  range. 

§ 
Modelling,  for  example,  combat  between  two  weapon  systems  with  different  maximum 

effective  ranges. 
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Section  6)  and  illustrated  these  theoretical  results  by  considering  some  numerical 
examples  (see  Section  7). 

Curves  of  the  modified  parity-condition  parameter  Z   plotted  against  the 
modified  offset  parameter  y      such  as  those  shown  in  Fig.  A  allow  one  to  parametrically 
analyze  "modern"  combat  modelled  with  the  general  power  attrition-rate  coefficients 
(9).   For  example,  we  can  now  parametrically  determine  whether  the  defender  will 
be  overrun  in  Bonder's  [1,3,8]  constant-speed-attack  model  (10)  with  attrition-rate 
coefficients  (11)  without  having  to  compute  the  entire  force-level  trajectories.   We 
illustrated  this  analysis  capability  with  some  numerical  examples,  which  showed  that 
the  defender's  annihilation  (i.e.  saturation  of  his  defensive  position  with  offensive 
fire)  depended  on  the  initial  force  ratio  (of  defender  to  attacker)  being  below  a 
certain  threshold  value.   Our  new  results  allow  one  to  not  only  easily  determine  such 
force-ratio  thresholds  of  survivability  but  also  study  their  dependence  on  weapon- 
system-capability  parameters. 

Our  new  results  let  us  conveniently  obtain  much  valuable  information  about 
the  model  (7).    Previously  one  was  limited  to  only  computing  force-level  trajectories, 
but  now  we  can  predict  battle  outcome  (in  particular,  force  annihilation)  without 

explicitly  computing  force-level  trajectories.   Moreover,  these  new  results  facilitate 

§ 
parametric  analysis   of  such  combat  situations.   In  particular,  Theorems  1  and  2 

explicitly  exhibit  a  tradeoff  between  quality  (as  quantified  by  the  relative-f ire- 

* 
effectiveness  parameter  X        and  the  parity-condition  parameter  Q  )   and  quantity 

(as  quantified  by  the  initial  force  ratio  xn/yn)  of  two  weapon  systems  in  combat 

against  each  other.   In  other  words,  one  can  use  an  expression  like  (21)  to  develop 


Varying,  for  example,  the  maximum  effective  range  of  the  defender's  weapons. 

The  classic  ordinary  differential  equation  theories  (see,  for  example,  Ince  [15]) 
were  inadequate  to  answer  many  important  questions  (for  example,  "Who  will  win? 
Be  annihilated?")  about  such  combat  models. 

§ 
S.  Bonder  [5]  has  suggested  that  an  increased  emphasis  be  placed  on  parametric 

analyses  in  systems  analysis  studies  (see  pp.  21-22  of  [5]). 
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quantitative  insights  into  how  the  quality  of  a  weapon  system  may  be  substituted 

for  sheer  numbers.   Moreover,  an  unanswered  theoretical  question  is  to  determine 

how  the  parity-condition  parameter  Q   depends  on  the  combat-intensity  parameter 

A,  and  the  relative-fire-effectiveness  parameter  X   .      Finally,  our  results  here 
I  k 

are  signposts  as  to  the  difficulty  of  analytically  extracting  information  (particu- 
larly parametric  information  without  excessive  computations)  from  variable-coef- 
ficient Lanchester-type  models  such  as  (7) . 
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